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Oh, Abstract 
(D . 

' We discuss the Next-to-Minimal Supersymmetric Standard Model (NMSSM) with a Peccei- 

^ ' Quinn (PQ) U{1) symmetry. When this symmetry is dynamically broken by the Higgs niech- 

' anism, the resulting pseudo-Nambu-Goldstone boson takes the role of an axion. Although 

' much of the allowed parameter space for low values of the PQ scale has been ruled out, many 



scenarios with a PQ scale > 10^ GeV remain untested, allowing the NMSSM PQ axion to 
provide a solution to the strong CP problem and be a good dark matter candidate. Un- 
fortunately the new particle states are so decoupled that they would not be observable at 
future colliders, and the NMSSM would appear indistinguishable from the minimal model. 
However, we show that in order to maintain vacuum stability, such a model requires that the 
heavy Higgs boson states have masses that lie close to approximately /xtan/3. Therefore, a 
measurement of the Heavy Higgs boson masses at the LHC would allow one to either rule out 
the NMSSM PQ axion, or provide tantalizing circumstantial evidence for its existence. 



Introduction: The Strong CP Problem and the Axion 



For some some time after its formulation, one of the principle strengths of Quantum Chro- 
modynamics (QCD) was thought to be its automatic conservation of parity (P) and charge- 
conjugation-parity (CP) symmetries. The only renormalizable P and CP violating term that 
may be added to the QCD Lagrange density is the "^-term" , 

A = e^^^F'^'^'^F;,, (1) 

where F^^°' is the gluon field strength and F^^ = ^e^i/po--^'"^" is it dual; ^eflf is the effective 
0-parameter after diagonlization of the quark mass matrix, i.e. ^eff = + arg det Mq. It is 
straightforward to show that this term is a total derivative allowing its integral over space-time 
to be written as a boundary term at infinity. Therefore, it was thought, its integral will vanishes 
in the vacuum, and the 0-term may be safely ignored. 

However, it was soon realized that such a term could not be ignored if the vacuum has non- 
trivial topological structure [1-3]. Indeed, even if set to zero by hand in the QCD Lagrange 
density, it will be regenerated when contributions from instanton solutions are included in the 
path integral. Its space-time integral does not necessarily vanish but is proportional to the 
winding number (Pontryagin index) of the field configuration. The 0-term will then contribute 
intrinsically non-perturbative CP violation, i.e. its effects will be invisible to perturbation theory. 
Since no CP violation has been observed in QCD, ^eflf must be very small. 

This can be quantified by examining the electric dipole moment of the neutron, d„: 
the CP violation induced by the 0-term leads to a neutron electric dipole moment of or- 
der ~ |0efr|lO~^^e cm which must be compared to the current experimental limit 
dn < 0.63 X 10~^^e cm [3]. Therefore |0eff| ^ 10~^, naturally leading to the question: why is 
CP violation in QCD so small? This is known as the "strong CP problem". 

The axion provides a very natural solution to the strong CP problem. It was realized that 
the 0-term could be absorbed by making a redefinition (an axial rotation) of the quark fields 
If the quarks have zero mass the Lagrange density will be unchanged except for the removal of 
the ^-term, and theories with differing values of ^^fr all represent the same physics. In essence, 
the ^-term can be rotated away using the global U(l) axial symmetry of the model. However, if 
the quarks have non-zero mass then this rotation will introduce complex phases into the quark 
mass matrix and the theory will still be CP-violating. 

Peccei and Quinn pj pointed out that if a new global axial symmetry, a Peccei-Quinn (PQ) 
symmetry, is introduced then it could be used to remove the ^-term instead. When this PQ 
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symmetry is dynamically broken by the vacuum structure it will result in a pseudo-Nambu- 
Goldstone boson known as the axion It is only a "pseudo"-Nambu-Goldstone boson because 
the PQ symmetry is not exact — it is explicitly broken by the triangle anomaly providing a 
non-perturbative axion-gluon coupling. This axion-gluon coupling has two effects. Firstly, it will 
provide a non-zero axion mass due to mixing with the pion, which is approximately given by 

' 'lOiOGeV" 



= [1 + 0(m„,,/m,)] « 0.6 x IQ-^eV 



f , (2) 

Ja 

where niu, rrid and mg and the up, down and strange quark masses respectively, /tt and m-,^ are 
the pion decay constant and the pion mass, and 0a is the axion field. Secondly, the axion-gluon 
coupling introduces an effective term into the Lagrange density of the same form as the 0-term, 
Eq.(^, so that the CP-violating terms become 

Ceof, = {e,^-^^^F^'"'F;,, (3) 

where fa is the axion decay constant. However, the potential for (pa is also a function of 
(^cff — (pa/ fa) and so the axion field relaxes to a vacuum-expectation- value (VEV) given by 
{(pa) = fadeS- The ^-tcrm is canceled and the strong CP problem is solved. 

The experimental bounds on the existence of the axion are already rather strict ;8;. The 
non-observation of an axion in collider experiments and rare decays (e.g. quarkonium decays) 
rules out models where the PQ scale (fa) is of the order of the electroweak scale. However, these 
bounds can always be avoided by increasing the PQ scale [HI EI| , or equivalently reducing the 
axion mass, thereby reducing the axion's couplings to known particles. 

In order to constrain this "invisible axion" one must consider astrophysical constraints 
Since a low mass axion is expected to be emitted during star cooling, fa may be constrained 
by insisting that the axion does not significantly alter the observed stellar evolution. Stars in 
globular clusters are the most sensitive to these effects [12-13]. Additionally, the neutrino signal 
from SN 1987A indicates that it is cooled mainly by neutrino emission rather than by emission 
of an "invisible axion" ^1]. Together these observations place a limit of roughly fa ^ 10^ GeV 
(translating via Eq.© to Ma < 0.01 eV). 

Intriguingly, at scales just above this limit the axion is seen to be a good dark matter can- 
didate. Indeed, it was shown in Ref.|15j that in the standard thermal scenario, and many in- 
flationary models, the dark matter axion 's PQ scale is predicted to be 3 x 10^° GeV. If 
the PQ scale becomes too much larger the axion contribution to dark matter may become too 
great, thereby over-closing the universe and thus providing an upper limit on fa- However, this 
upper bound is very model dependent. We will see later that the main results of this letter do 
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not depend on the fine details of the axion mass limits, but only that the PQ scale be very large. 

In this letter, we will briefly describe the PQ symmetric Next-to-Minimal Supersymmetric 
Standard Model (NMSSM), which is the minimal supersymmetric extension of Standard Model 
that can provide an axion. We will examine the Higgs boson mass spectrum of the model and 
see that the lightest pseudoscalar Higgs boson is the "invisible axion" , and will subsequently be 
unobservable at colliders for the foreseeable future. However, we will show that in order to keep 
the mass-squared of the lightest scalar Higgs boson positive, one must constrain the heavy Higgs 
bosons to lie in a very specific mass window. We will provide one-loop expressions for this mass 
window in a very good approximation. Therefore, this model provides a prediction for the heavy 
Higgs boson masses which may be confirmed or ruled out at the next generation of colliders. 

The PQ Symmetric NMSSM 

One model that provides an axion is the PQ symmetric NMSSM [15-18]; this has the same 
field content as the Minimal Supersymmetric Standard Model (MSSM) except for the inclusion 
of an extra Higgs singlet superfield S. Its superpotential is given, in an obvious notation, by 

W = u' h^QHu - hdQHd - h^LHd + \S{HuHd) ■ (4) 

The usual Higgs-higgsino mass term fiHuHd seen in the MSSM has been replaced by the term 
\S{HuHd) coupling the new singlet Higgs field, S, to the Higgs doublets, and H^^ where A 
is a dimensionless parameter. The Higgs-higgsino mass term will be recovered when the scalar 
component, 5, of the new singlet superfield gains a VEV of (5) = /i/A. 

In the MSSM, the dimensionful parameter, /x, is constrained to be of the order of the elec- 
troweak scale in order to give the correct pattern of electroweak symmetry breaking, even al- 
though it has no a priori relation to the electroweak scale. The question of why two seemingly 
unrelated scales should be the same is known as the "/i-problem" jTlj . The original formulation 
of the NMSSM was intended to answer this question by dynamically linking the scale to a VEV 
of a Higgs field, S*, and thereby to the electroweak scale. 

The superpotential, Eq.Q), has no dimensionful couplings and exhibits a U(l) PQ symmetry, 
which will be carried over into the Lagrange density. In the MSSM this PQ symmetry is explicitly 
broken by the Higgs-higgsino mass term ^HuHd', in the PQ symmetric NMSSM the PQ symmetry 
is only dynamically broken when S gains a non-zero VEV, giving rise to a near massless pseudo- 
Nambu-Goldstone boson — the axion^. Therefore the PQ symmetric NMSSM is the minimal 
^The axion is only a "pseudo"-Nambu-Goldstone boson since the PQ symmetry is explicitly broken by the 
triangle anomaly, giving it a small mass, Eq.lj^J. 
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supersymmetric extension of the Standard Model that can provide an axion. In fact, it is a 
supersymmetric version of the DFSZ axion model |lUj . 

The axion constraints mentioned in the introduction must also be applied here and so models 
with (S) of the order of the electroweak scale are ruled out. In the more usual formulation of 
the NMSSM this is avoided by adding a term ^kS^ to the superpotential; this explicitly breaks 
the PQ symmetry, giving the 'axion' a mass and avoiding the constraints. Here, in order to 
preserve a near massless axion, we insist that (5) > 10^ GeV. Therefore, the PQ symmetric 
NMSSM no longer links {S) to the electroweak scale and cannot be considered as a solution to 
the /i-problem. Since /i must remain of order the electroweak scale, A = /u/(5) must be very 
small and the /i-problem is re-expressed as: why is A so small? We will not attempt to answer 
this question here. 

The axion within the context of the NMSSM has also been discussed in Ref.jl9j. In that study, 
the term ^kS^ was included in the superpotential, explicitly breaking the PQ symmetry, but it 
was pointed out that in the limit where the soft supersymmetry breaking parameters associated 
with A and k vanish, the model will contain an additional approximate U{1)r symmetry. This 
symmetry is dynamically broken by the vacuum, giving rise to an 'R-axion'. Unfortunately the 
mass of this R-axion becomes rather large, forbidding its use in solving the strong CP problem, 
but nevertheless the model has interesting phenomenological consequences. 

The superpotential, Eqn.ljlJ), leads to the tree-level Higgs potential [T7| : 

V = VF + VD + V,olt, (5) 

with 

Vf = \XS\^{\H^\^ + \Haf) + \XH^Hdf, (6) 

Vd = lfi\H£-\H^\r+'^g'\HiH,\^ (7) 

^soft = mljHu\^ + mljHd\^ + ml\S\^ + [XAxSHuHd + h.c.], (8) 

where g = \/ g"^ + g''^ with g and g' being the gauge couplings of SU{2)l and U{1) interactions 
respectively, and we have adopted the notation HuHd = eap{Hu)°^{HdY = H:^H^ — H^H^. The 
first two terms, Vp and Vd, are the F and D terms derived from the superpotential in the usual 
way, while Vgoft contains the soft supersymmetry-breaking parameters A\, mn^, fnHi and ms- 
The vacuum of the model may be rendered neutral by a suitable application of a SU{2)l x 
U{1)y gauge transformation, and rendered real by exploiting the PQ symmetry. The vacuum is 
then given by 
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with Vs, Vu, and Vd real and positive. The requirement for this vacuum to be a local minimum 
provides three relations, linking the three soft mass parameters to the three VEVs of the Higgs 
fields: 

2 ^-2/ 2 2\ 1 \2 2 ) 1 /I \ 1 \2 2 

= gf K - Vd) - 2^ ^« + ~ 2 " ^ ^ 

2 -'--2/2 2\ ^ \2 2 , ^ A \ Vd 1,22 m\ 

= g9 [Vd - vj - 2^ ^rf + 71"^^^^^" ~ 2 " ^ ^ 

2 1x22, ^ \ A VuVd r,N 

= --A. +^AA--; (12) 



as usual, we have written v = \J Vu + v^^. 

The extra singlet fields mix with the Higgs doublet fields, increasing the rank of the scalar 
and pseudoscalar mass-squared mixing matrices by one each. After an initial rotation of the 
Higgs doublet fields by an angle /3, defined as usual via tan/3 = Vu/vd and outlined in detail in 
Ref . |i20] ■ the 2x2 pseudoscalar mass matrix is given by 

o / 1 isin2/3cot/3s \ , , 

V i sin 2/3 cot /3, i sin^ 2/3cot2 /?, j 

In analogy to tan/3, we have also defined tan/?^ = Vg/v; due to the requirement that 
{S) > 10^ GeV, tan /3s will be very large, and therefore cot /3s very small. In the above, we 
have defined the upper-left entry of the pseudoscalar mass-squared mixing matrix to be M\. 
This new mass parameter replaces the soft supersymmetry-breaking parameter Ax and becomes 
the mass of the MSSM pseudoscalar Higgs boson as the MSSM limit is approached, i.e. cot /3s ^ 
with fixed. This treatment allows higher order loop corrections to be absorbed directly into 
the definition of Ma- Including one- loop top/stop corrections, it is related to A\ by 

Ml = -^A, - ^AF{ml,ml) (14) 

where ht = ^/2mt / {v svn (3) is the top-quark Yukawa coupling and At is its associated soft 
supersymmetry-breaking mass parameter. The function F is given by 

nml.ml) = K K/^') - "^1 H/Q')] - 1 (15) 

and mt, rn^_^, m^^ are the top and stop masses, with Q the renormalization scale. 



This pseudoscalar mass-squared matrix is easily diagonalized, revealing two mass eigenstates, 
which will be denoted Ai and A2 with the label assigned in order of increasing mass. The first 
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of these, Ai, is the massless Nambu-Goldstone boson associated with the dynamical breaking of 
the PQ symmetry — the axion. The PQ symmetry ensures that it will be massless even after the 
inclusion of loop corrections; it only gains a very small mass via non-perturbative mixing with 
the pion, as described earlier. The heavier mass eigenstate, A2, has mass 

Ml, = + ^ sin2 2/3 cot^ /?,). (16) 

Since cot Ps ^ 10~'^ GeV the heavy pseudoscalar Higgs boson reproduces the mass of the MSSM 
pseudoscalar with a deviation less than one part in lO^"^. 



Similarly, the symmetric 3x3 scalar Higgs mass-squared matrix is 



+ A, 



(17) 



where the entries of the tree- level contribution, Mq, can be written as 



[Ml, 



12 



[^0]l3 
[MI22 

[MI23 



[Ml 



33 



Ml + (Ml; - cot^ (3s) sin^ 2(3 
-^(M|-/i2 cot2/3,)sin4/? 
1 9 

--Mjsin4^ cot/?s 

M| cos^ 2(3 + fi^ cot^ Ps sin^ 2P 

^(V - M^sin^ 2P) cot /3s 

^M^sin^ 2Pcot'^ps 



(18) 
(19) 

(20) 
(21) 
(22) 

(23) 



A denotes higher order corrections to the scalar Higgs mass matrix |221 123j . Including one-loop 
top/stop corrections these are given by [23] 



A 
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87r2 * 
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87r2 * 
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4 log 
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4 



8^Ki (m| , m| ) + 802^2 {ml , m? 
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'13 



S/3 



(24) 



,(25) 



levr^ 



o2 



-4:mfa—K2iml,ml) 



V2fi cot Ps 



3/1? 2 



2 2 
4 



8a'^b'^iyKi(mf ,m? ) -h 8ab^K2(mj ,mj 

*/3 



(26) 
(27) 
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A 



23 



3hl 
167r2 



3h 



—Am^acpK2 (in-^ , rii-^ ) \/2/x cot /3s 



A33 = -^m^Va c^cot^/3.^lHv"^|), 



(28) 
(29) 



where = sin/3, 0^3 = cos (3, a = (— /Ltc^ + Ats^)/\/2, and b = {jis 13+ AtCfj) / \/2, and the functions 
Ki and are 



2 \2 



K2(m|,m|) = (K(m?,m|) + l)/(m|+m|), 



with 



1 



(m^-^ , m^-^ ) = F{m^^ , m^-^ ) - - log 



2 2 



(30) 
(31) 

(32) 



Closed form expressions for the scalar Higgs boson mass eigenvalues can be obtained by diag- 
onalizing M^. However, these results are rather lengthy and unilluminating, and will not be 
reproduced here. 

Fortunately, these exact expressions are not needed due to the very small size of 
cot /3s ~ O(10~^). Notice that the mass-squared matrix takes the form 

All A12 Cicot/3s ^ 

M^= A12 A22 C2Cot/3s • (33) 

Ci cot Ps C2 cot /3s B cot^ /3s J 

This is true not only at tree- level but also when higher orders are included. We may reduce this 
matrix to block diagonal form by applying a unitary transformation defined by the 3x3 matrix 



1 - ^ cot2 /3srrt - cot /3sr 
cot/3srt 1- icot^^sF'^r 



+ 0(cot4/3s). 



with 



r = (Ci^22 - C2^i2, -C1A12 + CiAn) /det A. 
Applying this transformation gives the simple form 



(34) 



(35) 



( ^ii + Ciricot2/3s ^i2 + ^(Cir2 + C2ri)cot2/3s o ^ 

^12 + 5(Cir2 + C2ri)cot2/3s ^22+C2r2COt2/3s (36) 

\ (B - C2r2 - CiFi) COt2 /3s y 

+C'(cot3/3s). 
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The upper-left block consists of the usual MSSM scalar Higgs boson mass-squared matrix 
(to any desired number of loops) plus corrections of order cot^/3s. Consequently, the two heav- 
ier states, H2 and H^,, are rather uninteresting; the MSSM scalar Higgs masses, like a heavy 
pseudoscalar, will be recovered with corrections of only one part in 10^^, which is neither ex- 
perimentally observable, nor theoretically reliable since unincluded higher order corrections will 
present much larger deviations. This was to be expected since our NMSSM parameter choice is 
approaching the MSSM limit. 

A prediction for Ma 

The lightest Higgs boson, whose mass-squared is given by the lower-right entry, is rather more 
interesting. Its mass is suppressed by cot (3s , making it effectively massless at current collider 
energies, but its couplings to known particles, which mainly arise from the mixing with the other 
scalar Higgs bosons'^, will also be tiny. Subsequently, this state would be unobservable at high 
energy colliders for the foreseeable future, and the low energy phenomenology would appear 
indistinguishable from the MSSM. 

However, the expression for the lightest scalar mass shows interesting structure. Inserting 
the tree- level values into the lower-right entry of Eq. (|36|) gives the tree level mass-squared 

Mjj^ = fi^ tan^ 2(3 cot'^ (3si-{x^ + y^){x^ -if +y^cos^ 2(3)/ {xyf (37) 

where x = sin 2/3/(2//) and y = Mz sm2(3/{2fj,). 

This mass-squared must be positive in order to have a physically acceptable theory. If it is 
negative, the Higgs potential will be unbounded from below and the vacuum unstable. However, 
only the last term in the brackets of Ea. (|37() is positive; Mfj^ will become negative for both high 
and low values of Ma, and a stable vacuum will be achieved only for a small range around x ~ 1. 
This is also true when loop corrections are included, as shown in Fig.(^ left). 

To demonstrate this we examined 10^ different scenarios, with Ma and tan (3 chosen randomly 
between to 6 TeV and 3 to 30 respectively, we calculated the one-loop mass spectrum and, 
for every scenario with a stable vacuum, plotted a single point on the M^-tan /3 plane of Fig.(^ 
right). We discarded scenarios with unstable vacua. It is immediately evident that the physically 
acceptable scenarios all lie within a small band around Ma ~ 2/x/ sin 2/3 ~ /itan/3. 

Therefore the PQ symmetric NMSSM with a large expectation value of the new singlet field 
makes a prediction for the masses of the heavy Higgs bosons. This prediction is potentially falsifi- 
able, or verifiable, at the next generation of colliders. Furthermore, as long as cot Ps is small, the 
^The Lagrangian of the model also contains new direct couplings of the new singlet state to known particles 
but these are also suppressed by at least one order of cot /3s. 
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Figure 1: Left: The dependence of the lightest scalar Higgs mass (normalized by tanPs) on Ma, 
for tan/? = 10 and fi = 200 GeV . Beyond the points where the curve meets the axis the mass- 
squared becomes negative and the vacuum unstable. Right: The distribution of scenarios with 
physically acceptable vacua, with Ma chosen randomly between and 6 TeV and tan/3 chosen 
randomly between 3 and 30. The vacuum structure constrains the value of Ma to lie close to 
approximately /xtan/3. The blow-up allows individual scenario points to be seen. 



positivity or negativity of is independent of cot /3s, and consequently the prediction of the 
heavy Higgs boson masses is also independent of the value of cot Ps- Therefore, if after measuring 
fi and tan /3 at a future cohider, the heavy pseudoscalar mass is not found to he close to /i tan /3 
then this model is ruled out for all large values of the singlet expectation value. Alternatively, 
if the mass prediction were found to hold, it would provide very tantalizing, albeit indirect, ev- 
idence for the PQ symmetric NMSSM as a solution to the strong CP problem and for the PQ 
axion as a source of dark matter. 

In order to compare the values of /i and tan/3 with Ma at the next generation of colliders, 
the vacuum stability bounds must be made more precise. Since Mfj_^ = with M^^ given by 
Ea. ()37() is only a cubic in x^, it can be solved to give closed form analytical expressions for the 
tree-level boundary. Throwing the third non-physical solution away, we find 

^Lx/min = 1 - + - COS 7 ± V3smj) + A±, (38) 

where 

_ 1 _i / ^4(1 + 7/2)6 _ (2(1 +^2)3 _ 27^2 cos2 2^ 
^ ~ 3 2(l + y2)3_27y2 cos2 2/3 

and A-t represents the higher order corrections. 
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Since the one- loop top/stop contributions to A±, Eas. H24ll29|) . are independent of Ma, 
Mfj_^ = remains cubic in when these corrections are included and we can still find a closed 
form solution for the limits. However, these expressions are long and complicated, and once again 
such complexity is not needed here. Instead we expand the one-loop corrections as a series in 
the small parameter y and discard terms of 0{y^). This gives 

1/8 '-'^ 

^± = 7^ — A22 + [32A23 + 8si^A22 - 16S2/3C2,3Ai2] 

=Fy^^^ [8S2/3C2/3A1I - 852/3(1 + slp)c2pAl2 + 4/3(3 + S2/3)A22 

•S2/3C2/3 

-32s2/3C2^Ai3 + 16S2/3C2;3A23 + 32C2^A33] 

+16y2 [slpclp^ii + s^^C2;3Ai2 - 2S2/3C2/3A13] ^ + 0(y3), (40) 

where Ajj are given by Eas. piH^ . 

This approximation is rather good. The non-observation of supersymmetry to date requires 
that tan/? > 3 and /U > 80 GeV, giving y < 0.34. The discarded terms will therefore alter the 
one- loop corrections by at most a few percent. For more typical MSSM parameter choices, y will 
be even smaller; e.g. for the Snowmass reference point SPS la j24j, tan/3 = 10 and /x « 350 GeV, 
giving y ~ 0.026 

A large tan (3 expansion of the tree-level result gives a very approximate, but rather useful, 
"rule of thumb" : 

iVl4 « //tan/3± Mz. (41) 



The coupling of the lightest scalar Higgs boson to electrons may also be restricted by as- 
trophysical data, allowing more stringent limits to be placed on the PQ scale. Just as for the 
axion, Hi will be produced during the cooling of globular-cluster stars if its mass is below about 
10 keV. The maximum value of the Hi mass seen in Fig.(^ left) is realized^ at x ~ 1; inserting 
this into Ea. (|S7j) gives 

M^f ^ « 11 sin 2/3 cot /?« , (42) 

so the limits from star cooling cannot be avoided if {S) > 2/_isin2/3 x 10^ > 10^° GeV, where for 

the last inequality we have made the reasonable assumption that ^ ^ 1 TeV and tan /? > 3. 

^More accurately, making a series expansion in the small parameter y, the maximum (tree-level) value of Mhx 
is found at x = 1 + iy^ cos^ 2/3 + 0(1/"). 
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Above this scale one must respect the hmits on the couphng of the hghtest scalar Higgs boson 
to electrons [T^ . 

QH.e < 1.3 X 10-1^ (43) 

In the NMSSM it is easy to see that guie ~ nT-e/iS) and so this translates into a lower bound 
on the PQ scale. Combining this with the requirement that the Hi mass be less than 10 keV for 
this lower bound to apply, excludes the values 

2/xsin2/3 X 10^ < (5) < 4 X 10^°GeV. (44) 

Allowing the maximum and minimum values of /j, and tan f3 respectively, only a rather small range 
of {S) values is unequivocally ruled out. However, as /i and tan (3 are allowed to move toward 
less extreme values, the excluded range becomes larger and soon overlaps with that disallowed 
by emission of the axion from globular-cluster stars, i.e. (5) > 10^ GeV. 

Finally, since the model is supersymmetric, the extra neutral singlet superfield also contains 
a higgsino, which will be manifest as an extra neutralino — the lightest supersymmetric particle 
(LSP) of the model. Once again, the large value of the PQ scale leads to it having a very small 
mass and being almost totally decoupled from the other particles. To a good approximation, its 
mass IS given by Mlsp ^ fj. cot^ (3s, which, for /x « 10^ GeV and {S) « 10^^ GeV translates to^ 
-^LSP ~ 3 X 10^^ eV. In contrast to the scalar and pseudoscalar Higgs bosons, R-parity conserva- 
tion prevents the LSP being emitted during star cooling, so it provides no further astrophysical 
limits. 



Summary & Conclusions 

In this letter, we have discussed the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) with an explicit Peccei-Quinn (PQ) symmetry. This model is the minimal super- 
symmetric extension of the Standard Model that can provide an axion. This axion is a pseudo- 
Nambu-Goldstone boson associated with the dynamical breaking of the PQ symmetry, and is 
manifest in this model as the lightest pseudoscalar Higgs boson; it can be used to solve the strong 
CP problem of QCD and is a dark matter candidate. The stellar evolution of globular cluster 
stars and the neutrino signal from SN 1987A provide a lower bound on the PQ breaking scale 
> 10^ GeV. 

We have shown that in this limit simple expressions for the NMSSM Higgs boson masses 
can be obtained. The heavy and intermediate mass Higgs bosons have masses and couplings 
^It is intriguing to note that this mass Ues not too far from the expected neutrino mass scale. 
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indistinguishable from those of the corresponding MSSM. The hghtest scalar and pseudoscalar 
(the axion) decouple from the other particles and will be invisible to future collider searches. 

However, we have demonstrated that in order that the theory have a stable vacuum, i.e. in 
this case that the lightest scalar mass-squared be positive, the heavy mass scale Ma must lie 
within approximately ^tan/3 it Mz- We have presented analytic expressions for these limits on 
Ma to one- loop top/stop accuracy. 

If, at a future collider. Ma were found to be outside this range, then the PQ symmetric 
NMSSM would be ruled out for all values of the PQ scale. This is not an unreasonable event; 
the restriction on Ma is unlikely to occur by chance without some other organizing principle. For 
example, all of the Snowmass MSSM reference points 24 , which are considered a representative 
sample of MSSM scenarios, fail this criterion. It is important to stress that only the axion 
associated with this particular model would be ruled out; an axion could still be present via 
some other mechanism, and axion search experiments, such as CAST |25j . the U.S. Axion Search 
(Livermore) j26j and the Kyoto search experiment CARRACK would still be very important. 

On the other hand, if the heavy Higgs boson mass scale were seen to obey the bound given 
by Ea. (|38|) we would have very exciting circumstantial evidence for the existence of an NMSSM 
axion. Then the role of the axion search experiments would become even more crucial. 
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